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Solving the Gross-Pitaevskii-type equation it is shown that the magneto-oscillations, observed 
in the superconducting state of a few underdoped cuprates, originate in the quantum interference 
of the vortex lattice with nanoscale crystal lattice modulations of the order parameter as revealed 
by Scanning Tunnelling Microscopy (STM). The commensuration of the vortex lattice and crystal 
lattice have 1/ B^^^ periodicity, rather than 1/B periodicity of conventional normal state magneto- 
oscillations. Experimental conditions allowing for a resolution of two different types of quantum 
magneto-oscillations are outlined. 

PACS numbers: 74.20.-z,74.72.-h,71.38.-k 
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Until recently no convincing signatures of quantum 
magneto-oscillations have been found in the normal 
state of cuprate superconductors despite significant ex- 
perimental efforts. There are no normal state oscilla- 
tions even in high quality single crystals of overdoped 
cuprates like Tl2Ba2Cu06, where conditions for de Haas- 
van Alphen (dHvA) and Shubnikov-de Haas (SdH) os- 
cillations seem to be perfectly satisfied [i| and a large 
Fermi surface is identified in the angle-resolved pho- 
toemission spectra (ARPES) [2|- The recent observa- 
tions of magneto-oscillations in kinetic 0, lj| and mag- 
netic '5^ response functions of underdoped YBa2Cu306.5 
and YBa2Cu408 are perhaps even more striking since 
many probes of underdoped cuprates including ARPES 
p] clearly point to a non Fermi-liquid normal state. Their 
description in the framework of the standard theory for 
a metal [Tf has led to a very small Fermi-surface area 
of a few percent of the first Brillouin zone [3, 0, Q , and 
to a low Fermi energy of only about the room tempera- 
ture Q- Clearly such oscillations are incompatible with 
the first-principle (LDA) band structures of cuprates, but 
might be compatible with a non-adiabatic polaronic nor- 
mal state of charge-transfer Mott insulators Q. Never- 
theless their observation in the suverconducting (vortex) 
state well below the Hc2iT)-\me [3] raises a doubt con- 
cerning their normal state origin. 

Here I propose an alternative explanation of the 
magneto-oscillations 0, 0, [^ as emerging from the quan- 
tum interference of the vortex lattice and the checker- 
board or lattice modulations of the order parameter ob- 
served by STM with atomic resolution [9|. The checker- 
board effectively pins the vortex lattice, when the pe- 
riod of the latter A = {irh/eBY/'^ is commensurate with 
the period of the checkerboard lattice, a. The condition 
A = Na , where TV is a large integer, yields 1/ B^'^ peri- 
odicity of the response functions, rather than 1/B period- 
icity of conventional normal state magneto-oscillations. 

To illustrate the point one can apply the Gross- 
Pitaevskii (GP)-type equation for the superconducting 
order parameter V'(r), generalized by us [lO| for a charged 
Bose liquid (CBL), since many observations including a 
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FIG. 1: The order parameter profile /(p) = '4>{r)/ns of a 
single vortex in CBL [I^ (symbols). Inset: CBL vortex (a) 
[lOl . [lq | compared with the Abrikosov vortex (b) 14] (here 



small coherence length point to a possibility that under- 
doped cuprate superconductors may not be conventional 
Bardecn-Cooper-Schrieffer (BCS) superconductors, but 
rather derive from the Bose-Einstein condensation (BEG) 
of real-space pairs, such as mobile small bipolarons [11], 
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r')lV(r')l' ^(r)-0. 

(1) 

Here -E(K) is the center-of-mass pair dispersion and the 
Peierls substitution, K => —ihW + 2eA is applied with 
the vector potential A(r). The integro-differential equa- 
tion ^ is quite different from the Ginzburg-Landau [I4] 
and Gross-Pitaevskii [l^ equations, describing the order 
parameter in the BGS and neutral superfluids, respec- 
tively. In the continuum (effective mass) approximation, 
E(K) = h'^K'^ /2m**, with the long-range Coulomb re- 
pulsion between double charged bosons, ^(r) = ^(r) = 



4e^/eor, this equation describes a single vortex with a 
charged core, Fig.l, and the coherence length roughly 
the same as the screening radius, $, — {h/2^/'^m**u}pY/'^. 
Here ujp = (167rnse^/eom**)^/^ is the CBL plasma fre- 
quency, eo the static dielectric constant of the host lat- 
tice, m** the boson mass, and Us is the average con- 
densate density. The chemical potential is zero, /i = 0, 
if one takes into account the Coulomb interaction with 
a neutralizing homogeneous charge background, or de- 
fines the zero-momentum Fourier-component of V^(r) as 
zero. Each vortex carries one flux quantum, </'n_= 7r/i/e 
, but it has an unusual core. Fig. la, Ref. [lO], due 
to a local charge redistribution caused by the magnetic 
field, different from the conventional vortex [ij|. Fig. lb. 
Remarkably, the coherence length turns out very small, 
^ « 0.5nm with the material parameters typical for un- 
derdoped cuprates, m** — 10?Tie, Ug — lO^^crn^"^ and 
eo = 100. The coherence length ^ is so small at low 
temperatures, that the distance between two vortices re- 
mains large compared with the vortex size, A 3> ^, [l5| 
in any laboratory field reached so far 0, 0, Q- It al- 
lows us to write down the vortex-lattice order parameter, 
ii{v) = i/jviir), as 



ipvi{r) 



,1/2 



E< 



r,) 



(2) 



where (^(r) = 1 — f{p), and r^ = Aj? 
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0, ±1, ±2, ... (if, for simplicity, we take the square vortex 
lattice). The function (j){p) is linear well inside the core, 
(j){p) « 1 — 1.52p/^ (p ^ C)i and it has a small negative 
tail, (j){p) « —A^^/p"^ outside the core when p ^ ^, Fig.l 

In the continuum approximation with the Coulomb 
interaction alone the magnetization of CBL follows the 
standard logarithmic law , M{B) oc In 1/B without any 
oscillations since the magnetic field profile is the same as 
in the conventional vortex lattice [l4| . However, more of- 
ten than not the center-of-mass Bloch band of preformed 
pairs, £^(K), has its minima at some finite wave vectors 
K = G of their center-of-mass Brillouin zone [Il|, [iq . 
Near the minima the GP equation ([1]) is written as 
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with the solution ip{r) = V'G(r) = e V-ui(r), if the in- 
teraction is the long-range Coulomb one, V{r) = Vc(r). 
In particular, a nearest-neighbor (nn) approximation for 
the hopping of intersite bipolarons between oxygen p- 
orbitals on the Cu02 2D lattice yields four generate 
states "00 with G,; = {±27r/ao, ±27r/ao} , where oq is 
the lattice period (l6[. Their positions in the Brillouin 
zone move towards F point beyond the nn approximation. 
The true ground state is a superposition of four degen- 
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FIG. 2: The checkerboard d-wave order parameter of CBL 
[rn | on the square lattice in zero magnetic field (coordinates 
x,y are measured in units of a). 



erate states, respecting time-reversal and parity symme- 
tries [i3|, 

V'(r) = An]/'^ [cos(7rx/a) ± cos(7ry/a)] V^i(r)- (4) 

Here we use the reference frame with x and y axes along 
the nodal directions and a — 2^^'^ao. Two "plus/minus" 
coherent states, Eq.(l4|), are physically identical since they 
are related via a translation transformation, y => y + 
a. Normalizing the order parameter by its average value 
{tp{r)'^) — Us and using (f/A)^ ^ 1 as a small parameter 
yield the following "minus" state amplitude, A » 1 — 

(j)2{2TT/a)]dn^ii for the square vortex lattice |18j] with the 
reciprocal vectors g — {2TT/X){nx,ny}. Here Sn,R is the 
Kroneker symbol, R — A/a is the ratio of the vortex 
lattice period to the checkerboard period {n — 0, 1, 2, ...), 
N = BS/ (j)Q is the number of flux quanta in the area S 
of the sample, and (j)k{q) = i'^'n'/S) J^ dppJo{pq)(l)''{p) is 
the Fourier transform of fc's power of 0(p), where Jo{x) 
is the zero-order Besscl function. 

The order parameter ^/'(r), Eq.Q has the d-wave sym- 
metry changing sign in real space, when the lattice is 
rotated by 7r/2. This symmetry is due to the pair center- 
of-mass energy dispersion with the four minima at K 7^ 0, 
rather than due to a specific symmetry of the pairing 
potential. It also reveals itself as a checkerboard mod- 
ulation of the carrier density with two-dimensional pat- 
terns in zero magnetic field. Fig. 2, as predicted by us jl7[ 
prior to their observations 0. Solving the Bogoliubov- 
de Cennes equations with the order parameter, Eq.(|3]), 
yields the real-space checkerboard modulations of the 
single-particle density of states [13] , similar to those ob- 
served by STM in cuprate superconductors. 

Now we take into account that the interaction be- 
tween composed pairs includes a short-range repulsion 
along with the long-range Coulomb one, V{r) = 14 (r) -I- 



vS{r) [la|- At sufficiently low carrier density the short- 
range repulsion is a perturbation to the ground state, 
Eq-O, if the corresponding characteristic length, ^/^ = 
h/{2m**nsvy^'^ is large compared with the coherence 
length ^, related to the long-range Coulomb repulsion, 
C/i ^ C- The short-range repulsion constant v is roughly 
the pair bandwidth w of the order of 100 meV times 
the unit cell volume, v w waj] [Ig. Using this estimate 
one can readily show that the perturbation treatment 
of the short-range interaction is justified for any rele- 
vant density of pairs, if eg < 10"^. On the other hand, 
a strong short-range interaction could affect both the 
checkerboard and the vortex lattices, if (,h is compara- 
ble with ^ and a. 

Importantly the short-range repulsion energy of CBL, 
U = {v/2){ip{r)^), has a part, AU, oscillating with the 
magnetic field as 
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where Uq ~ 'vnl/2 is the 



hard-core energy of 

the amplitudes are 
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- 4502 (2i/2^/a) + 

-f 801 (lOi/^TT/a) - 

(57/2)02 (27r/a) - 
1201 (23/2^/a) + 



a homogeneous CBL, and 
proportional to the Fourier 
as Ai = 1501 (2i/27r/a) - 
2403 (2i/27r/a) - 604 (2i/2^/a) - 
1202 (l0^/^7r/a) + 803 (lO^/^Tr/a) 
A2 - -(23/2)01 (27r/a) + 
1603(27r/a) -f 404 (27r/a) - 
902 (23/27r/a)-603 (23/27r/a) -f 304 {2'^'^-K/a), and A3 = 
-01 (47r/a)-f (3/2)02 (4^/a)-03 (47r/a)-f (1/4)04 (Ma). 
Fluctuations of the pulsed magnetic field and unavoid- 
able disorder in cuprates induce some random distri- 
bution of the vortex-lattice period, A. Hence one has 
to average AU over R with the Gaussian distribution, 
G{R) = exp[-(i? — ^_)2/^2j^^^i/2 around an average R 
with the width 7 <C i?. Then using the Poisson summa- 
tion formula yields 
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cos(47rfci?). (6) 



The oscillating correction to the magnetic susceptibil- 
ity, Ax{B) — —d^d/dB^, is strongly enhanced due to 
high oscillating frequencies in Eq.Q. Since the super- 
fluid has no entropy we can use AU as the quantum 
correction to the thermodynamic potential fl even at fi- 
nite temperatures below Tc{B). Differentiating twice the 
first harmonic (fc = 1) of the first lesser damped term in 
Eq.® we obtain 



Ax(S) « xoe 



-S^Bo/16B 



{?^^\os{Bo/BY'\ (7) 
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FIG. 3: Quantum corrections to the vortex-lattice suscep- 
tibility versus 1/B, Eq.0 (solid line, Bo = 1-000 • 10^ 
Tesla, 5 — 0.06) compared with oscillating susceptibility of 
YBa2Cu306.5 (symbols) and with the conventional normal 
state oscillations (dashed line) [lyl at T = 0.4 K, upper panel. 
Middle panel: oscillating susceptibility of YBa2Cu4 08 (sym- 
bols [Sl) at T = Q.53K compared with the theory (solid line, 
Bo = 1-190 ■ 10** Tesla, S = 0.07), where 20 percent of the 
triangular vortex lattice susceptibility added to the square 
lattice one. Lower panel: quantum corrections to the current. 



where Xo = UoSA^e a /4n h is a temperature- p^portional to F^/B (Bo = 0.853 - 10« Tesla and S = 0.1 



dependent amplitude, proportional to the condensate 



solid line) compared with the oscillatory part of the Hall re- 
sistance in the mixed state of YBa2Cu306.5 (symbols [3[) at 
1.5 K. 



density squared, Bq = ■K^hjed} — ^Ti^h/ea^ is a char- 
acteristic magnetic field, which is approximately 1.1 • 10^ 
Tesla for oq ~ 0.38 nm, and 7 is replaced by 7 = 5R with 
the relative distribution width 5. Assuming that ^ 



> 



a, 



so that the amplitude Ai is roughly a? / S, the quantum 
correction Ax, Eq.l]?]), is of the order oi wx^ / B^ , where x 
is the density of holes per unit cell. It is smaller than the 
conventional normal state (de Haas-van Alphen) correc- 
tion, AxdHvA ^ fJ-/B^ [7], for a comparable Fermi-energy 
scale ^ — wx, since a; ^ 1 in the underdoped cuprates. 

Different from normal state dHvA oscillations, which 
are periodic versus 1/B, the vortex-lattice oscillations, 
Eq.lO are periodic versus 1/B^^^. They are quasi- 
periodic versus 1/B with a field-dependent frequency 
F — Bq^B/BoY^^ /27r, which is strongly reduced rel- 
ative to the conventional-metal frequency (w BnliTr) 
since i? <C -Bq, as observed in the experiments [3, 0j 0- 
The quantum correction to the susceptibility, Eq.© fits 
well the oscillations in YBa2Cu306.5 19], Fig. 3 (upper 
panel) . Importantly, if the vortex lattice has two domains 
with different coordination of vortices [l3|, then there 
are two resonating fields, -Bq of the square lattice and 
Bi = 2B0/3 of the triangular lattice, causing beats in the 
oscillations, as observed in Ref. [a| at low temperatures 
in YBa2Cu408(Fig.3, middle panel). A pinning force on 
the vortex lattice, Fp, due to the checkerboard modula- 
tions is proportional to dU /da. Hence the oscillating part 
of the Hall and longitudinal resistivity is proportional to 
Fp/B oc exp-52^o/16S(Bo/S)l/2sin(Bo/5)^/^ which 
fits the oscillatory part of the Hall resistance 3] rather 
well, Fig. 3 (lower panel). The oscillations amplitudes, 
proportional to n^ exp(— 5^_Bo/16B) decay with increas- 
ing temperature since the randomness of the vortex lat- 
tice, 5, increases, and the Bose-condensate evaporates. 

In summary, I propose that the magneto-oscillations 
in underdoped cuprate superconductors [3, 0, Q result 
from the quantum interference of the vortex lattice and 
the lattice modulations of the order parameter. Fig. 2, 



playing a role of the periodic pinning grid. The mag- 
netic length, A > 5 nm, remains larger than the zero- 
temperature in-plane coherence length, £, "^ 2 nm, mea- 
sured independently, in any field reached in Ref. 0, 0, • 
Hence the magneto-oscillations are observed in the vortex 
(mixed) state well below the upper critical field, rather 
than in the normal state, as also confirmed by the neg- 
ative sign of the Hall resistance J3J. It is well known, 
that in "YBCO" the Hall conductivities of vortexes and 
quasiparticles have opposite sign causing the sign change 
in the Hall effect in the mixed state [20]. Also there is 
a substantial magnetoresistance [J], which is a signature 
of the fiux fiow regime rather than of the normal state. 
Hence it would be rather implausible if such oscillations 
have a normal-stjite origin due to small electron Fermi 
surface pockets 19] with the characteristic wave-length 
of electrons larger than the widely accepted coherence 
length. In any case our expression [71 describes the oscil- 
lations as well as the standard Lifshitz-Kosevich formula 
of dHvA and SdH eflFects @, H H, [3. The difference 
of these two dependences could be resolved in ultrahigh 
magnetic fields as shown in Fig. 3, upper panel. 

While our theory utilizes GP-type equation for hard- 
core charged bosons [lOj . the quantum interference of 
vortex and crystal lattice modulations of the order pa- 
rameter is quite universal extending well beyond Eq.([T]) 
independent of a particular pairing mechanism. It can 
also take place in the standard BCS superconductivity at 
B < Hc2, but hardly be observed because of much lower 
value of Hc2 in conventional superconductors resulting in 
a very small damping factor, ex exp(— 5^i?o/16-B) <SC 1. 
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